Normal Mode Splitting and Antibunching in Stokes and Anti-Stokes Processes in 
Cavity Optomechanics: Radiation Pressure Induced Four- Wave Mixing Cavity 

Optomechanics 



O 



o 

m 

S 1 

> ■ 

c ■ 
ctf ■ 
=5 

cr 



> 

00 

ON 
O 

o 



X 



Sumei Huang and G. S. Agarwal 
Department of Physics, Oklahoma State University, Stillwater, Oklahoma 74078, USA 

(Dated: March 31, 2010) 

We study Stokes and anti-Stokes processes in cavity optomechanics in the regime of strong cou- 
pling. The Stokes and anti-Stokes signals exhibit prominently the normal-mode splitting. We 
report gain for the Stokes signal. We also report lifetime splitting when the pump power is less 
than the critical power for normal-mode splitting. The nonlinear Stokes processes provide a useful 
method for studying strong coupling regime of cavity optomechanics. We also investigate the corre- 
lations between the Stokes and anti-Stokes photons produced spontaneously by the optomechanical 
system. At zero temperature, our nanomechanical system leads to the correlations between the 
spontaneously generated photons exhibiting photon antibunching and those violating the Cauchy- 
Schwartz inequality. 

PACS numbers: 42.50.Wk, 42.65.Dr, 42.65.Ky 



I. INTRODUCTION 

The nonlinearities in a system can be studied using 
a number of optical methods. Among these, Stokes 
and anti-Stokes processes, and more generally four-wave- 
mixing processes are quite common tools used to under- 
stand the nonlinear nature of the system | lj . With this 
in view we study the stimulated Stokes and anti-Stokes 
processes in cavity optomechanics. As is well known, the 
nonlinearity in cavity optomechanics arises from the ra- 
diation pressure [2|-|7| on the moving mirror of the cavity. 
Thus, if the cavity is driven by a pump field of frequency 
oji and a Stokes field of frequency w s , then, due to ra- 
diation pressure, the output of the cavity would consist 
of fields at the applied frequencies uii and uj s and a gen- 
erated frequency 2w; — w s . While some previous works 
[H-[l(| have explored the Stokes and anti-Stokes processes 
in the context of parametric oscillation instability, here 
we show how such processes can be conveniently used to 
study the phenomena of normal-mode splitting [TlT - fl9j 
arising from the strong coupling between the cavity and 
the mechanical mirror. Further, the system can act as an 
amplifier for the Stokes field. Needless to say, we work 
in a domain which is below the instability threshold. 

Moreover, very interesting photon correlations between 
the Stokes and the anti-Stokes photons have been re- 
ported in atomic vapors under conditions of electromag- 
netically induced transparency [20]. Here we also dis- 
cuss the correlations between the photons created spon- 
taneously by the optomechanical system. The correla- 
tions are found to be nonclassical. 

The article is organized as follows. In Sec. II, we in- 
troduce the model, obtain the equation of motion for the 
oscillator and the cavity field, and solve it. In Sec. Ill, 
we calculate the output fields and thus obtain nonlin- 
ear susceptibilities for Stokes and anti-Stokes processes. 
In Sec. IV, we show that the Stokes field is amplified, 
and find very prominent normal-mode splittings in the 



output fields. Thus, stimulated Stokes and anti-Stokes 
processes provide us with a new tool for studying the 
strong coupling regime of optomechanics. We find that 
normal-mode splittings are especially pronounced in the 
two quadratures of the output fields. In Sec. V, we 
analyze the correlations between the spontaneously gen- 
erated photons in the four-wave-mixing processes in the 
optomechanical system. We show that such correlations 
are intrinsically quantum. 



II. MODEL: STIMULATED GENERATION OF 
STOKES AND ANTI-STOKES FIELDS 

We consider the system illustrated in Fig. [TJ in which 
the cavity consists of two mirrors separated from each 
other by a distance L. The front mirror is fixed and par- 
tially transmitting; the end mirror is movable and per- 
fectly reflecting. The cavity is driven by a pump field 
and a Stokes field obtained with lasers. Their frequencies 
are uji and uj s , respectively. We would assume that the 
Stokes field is much weaker than the pump field. A radi- 
ation pressure produced by momentum transfer will act 
on the movable mirror, which is modeled as a harmonic 
oscillator with mass m, frequency u m , and momentum 
decay rate 7 m . 

Considering a single-mode cavity w c , the Hamiltonian 
of the system in a frame rotating at the pump frequency 
lji is written as 



H = h(uj c - uJi)n c - huj m xn c Q + 



■(Q 2 



p2 



+ihei(c' ! - c) + ift^e-^^M - e*e l ^-" ! )*c]. 

(1) 

Here Q and P are the dimensionless operators represent- 
ing the oscillator's position and momentum, defined by 
Q = q^j2mu) rn /K and P — Py/2/ (mhuj m ) with [Q, P] = 
2i. In Eq. ([T]), the first term is the energy of the cav- 
ity field, n c = e'e is the number of the photons inside 



2 



<y s , a>i, 2a>i — u) s 



fixed 



cavity axis 
■* ► 



movable 



FIG. 1: Sketch of the studied system. A pump field with 
frequency uji and a Stokes field with frequency ui a enter the 
cavity through the partially transmitting mirror. The output 
fields c ou t have three components (uji , u> s , 2cji — u s ). No vac- 
uum fields are shown here because we are examining only the 
mean response. 
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Thus Eq. ((3} shows the cavity field (c)e~ lUHt has three 
components, oscillating at the input frequencies ui and 
lo s , and a new anti-Stokes frequency 2ui — co s . By sub- 
stituting Eq. (j3|) into Eq. ([2]), neglecting those terms 

and equating coefficients of 



containing ef, e* 2 , and \e s \ 2 
terms proportional to e~^" s 
tively, we find 
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the cavity and c and are the annihilation and cre- 
ation operators, respectively, for the cavity field satis- 
fying the commutation relation [c, c^] = I. The second 
term describes the nonlinear coupling of the movable mir- 
ror to the cavity field via radiation pressure, where the di- 
mensionless parameter \ = / u m ){w c / L)^J h/ (2rau m ) is 
the optomechanical coupling constant between the cavity 
field and the movable mirror. The third term corresponds 
to the energy of the movable mirror. The last two terms 
give the interactions of the cavity field with the pump 
field and the Stokes field, ei and e s are, respectively, the 
amplitudes of the pump field and the St okes field in side 
the cavity. They are defined by e; = ^/2np/ (tiuii) and 
|e s | = y /, 2Kp s /(hio s ), respectively, where p is the pump 
power, p s is the power of the Stokes field, and k is the 
cavity decay rate due to the fixed mirror. 

Let (Q), (P), (c), and (c^) be the expectation values 
of the operators Q, P, c, and c\ respectively. The time 
evolution of these expectation values can be derived by 
using the Heisenberg equations of motion and adding the 
damping terms: 



(Q)=LU m (P), 

(P) = 2u> m x(n c ) - u) m (Q) - lm(P), 

(c) =-[* + i(u c -uji- u nX {Q))]{c) +ei+ ese-^'-^, 

(ct) = [k i{uj c -ui- u mX {Q)W) +ei+ eje^'-*)*. 

(2) 

The derivation of Eq. (|2|) uses the well-known mean-field 
assumption (Qc) — (Q)(c). As the field e s at the Stokes 
frequency u s is much weaker than the pump field we 
derive the steady-state solution of Eq. ^ to first order 



Qo = 2 X |c | 2 , 
Po = 0, 



co = 
c+ = 



k + iA ' 
1 



-{[k - i(A + u s - lui)] 



(4) 



(5) 



d(u s - LUl) 

x [{u s - ui) 2 ~u 2 m + ij m (u s - uii) 
-2iu 3 mX 2 \c \ 2 }, 

d*(u s - lji) ' 

where 

A = U c - UJi - UJmXQo, 

is the effective detuning, and where 
d{u s - ui) = Alo^x 2 A\c Q \ 2 + [(u a uj m ) 

x (u s -U)i- u m ) + i7m(w s - ui)] 

x [k + i ( A — u s + ui )] [k — i ( A + u s — cj; )] . 

For brevity we do not write explicit expressions for Q±, 
P±, etc. because we do not need these in the discussion 
that follows. 



III. THE OUTPUT FIELDS 

To investigate normal-mode splitting of the output 
fields, we need to find the expectation value of the output 
fields. Using input-output relation (2~lj (c ou t) +Ei/V2k + 
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e s e l (" s W! ^/a/2k = \/2/c(c), we can obtain the expecta- 
tion value of the output fields 

(c out ) = V2k~[c +e s e-^-^)* c+ +e%e i ^-^ t cJ\ 

(7) 

If we write (c ou t} as 

(cout) =c t + ese-'^-^Cs + ey^'- u ^c as , (8) 

where c; is the response at the pump frequency u>i, c s 
is the response at the Stokes frequency u„ and c as is 
the response at the four- wave-mixing frequency 2loi — uj s 
(anti-Stokcs frequency). Then we have 
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In the absence of the interaction between the cavity field 
and the movable mirror, one would expect the output 
fields to contain only two input components (wj and w s ); 
no four-wave-mixing component appears. We can get this 
result from Eq. © by setting x — 0, which gives 



ci 



2kei 



iA 
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12k 



(10) 



+ i(A - u s + ui{) 



'2k 



= 0, 



as expected. However, in the presence of the coupling 
with the oscillator (x ^ 0), [from Eq. ([5]), we have 
q 7^ 0, c s 7^ 0,c as 7^ 0], the output fields contain three 
components. The generated signal would exhibit reso- 
nances whenever w s = ± uj m . In addition, one would 
have the resonances produced by the cavity uj s — uji ± A. 
These resonances are, of course, expected. The normal- 
mode splitting would arise as a result of strong coupling 
X [l7Hl9l|. This is because the structure of the denom- 
inator in Eq. © depends on x- We next present the 
roots of Eq. ©. 

We use parameters which have been used in a recent 
experiment on the observation of the normal-mode split- 
ting in the fluctuation spectra [l7( | : the wavelength of 



the laser A = 2irc/uii — 1064 nm, L = 25 mm, m = 145 
ng, k — 2tt x 215 x 10 3 Hz, ui m = 2n x 947 x 10 3 
Hz, the mechanical quality factor Q' = u m /^ m — 6700, 
j m = 27rxl41 Hz, A = u) m . In this range of parameters, 
no parametric instabilities occur. 

Figure [5] shows the dependence of the real parts of the 
roots of d(uj s — in the domain Re(w s — uii) > on 
the pump power. Figure [3] shows the dependence of the 
imaginary parts of the roots of d(ui s — u>i) on the pump 
power. For a small value of the pump power, the real 
parts of the roots of d(uj s — u>i) have two equal values, so 
there is no splitting. However, there is lifetime splitting 
[22| as seen in the Fig. [3] If we increase the pump power 
to a certain value, the real parts of d(uj s — u>i) in the do- 
main Re(a; s — > begin to have two different values, 
and the difference between two real parts of the roots of 
d(uj s — u>i) in the domain Re(w s — loi) > is increased 
with increasing pump power. 
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FIG. 2: (Color online) The roots of d(o; s — oji) in the domain 
Re(cj s — uji) > as a function of the pump power p. 
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FIG. 3: (Color online) The imaginary parts of the roots of 
d(ui s — uji) as a function of the pump power p. 



IV. NORMAL-MODE SPLITTINGS IN THE 
OUTPUT FIELDS 

Before examining the normal-mode splitting in, say, 
the output anti-Stokes field, we examine Eq. Q in the 
traditional limit of nonlinear optics; that is, we find the 
form of anti-Stokes field to lowest order in 



4 



(« + iA) 2 [(uj s - lui + io m ){uj s —cji — uj m ) - i-/ m (uj s - W{)][k - i(A - uj s + wj)][k + i(A + cj s - wj)] ' 



which has resonances as discussed after Eq. (|10p and 
which is proportional to the pump power. 

We next discuss the normal-mode splitting in the gen- 
erated Stokes and anti-Stokes fields. It is useful to nor- 
malize all quantities to the input Stokes power p s . For 
simplicity, we assume e s to be real. For our plots we 
would give the output power at the Stokes frequency u s 
in terms of the input Stokes power 

G s = ^^l = \V2-.c s f, (12) 

and the two quadratures of the output fields at the Stokes 
frequency lu s in terms of the square root of the input 
Stokes power. Let us denote these normalized quadra- 
tures by v s and v s . These are defined as v s — \/2k c " 2 c? 

and v s — V2k Cs 2 ^ 3 ■ The quantity G s is the gain of the 
cavity optomechanical four- wave mixer. In Figs. [3]- |5J 
we have plotted v s , v s , and G s , respectively, versus the 
normalized frequency (lo s — u>i ) /uj m for different pump 
powers. The quadrature v s (v s ) exhibits absorptive (dis- 
persive) behavior. As is known, there is a phase change 
on reflection and that is why the quadrature v s shows 
absorptive behavior. The normal-mode splitting or the 
lifetime splittings are clearly seen depending on the in- 
put pump power in the quadratures v s and v a . The peak 
positions are in agreement with Fig. [2] for the case when 
the input pump power is such that normal-mode split- 
ting occurs. The behavior of net gain as a function of 
lj s is different due to the combination of absorptive and 
dispersive characteristics of the quadratures v s and v s . 
The gain shows normal-mode splitting for larger value 
of the pump power. Moreover, the maximum gain of the 
Stokes field is about 1.15. It should be borne in mind that 
the quadratures v s and v s can be obtained by homodyne 
measurement. 

Likewise, the output power at the anti-Stokes fre- 
quency 2ui — lj s in terms of the input Stokes power is 
given by 

npun-^M* _ ^ 2 _ (13) 

Ps 

For brevity, we only show in Fig. [7J the function G as 
against the normalized frequency (w s — loi) /ui m f° r sev_ 
eral values of the pump power. As can be seen in Fig. [71 
increasing the pump power can make the signal of four- 
wave mixing evolve from one peak to double peaks. It 
is also seen that the maximum value of G as is about 
0.15 and the output power at the anti-Stokes frequency 
(2u>i — lo s ) is much less than the output power of the 
Stokes field. However, for larger pump powers, the max- 
imum gain for Stokes and anti-Stokes fields are bigger. 
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FIG. 4: (Color online) The normalized quadrature v s plotted 
as a function of the normalized frequency (lj s — uii)/ui m for 
different pump power, p = 1 mW (solid curve), 6.9 mW 
(dotted curve), and 20 mW (dashed curve). 
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FIG. 5: (Color online) The normalized quadrature v s plotted 
as a function of the normalized frequency (u> s — oJi)/oJ m for 
different pump power, p — 1 mW (solid curve), 6.9 mW 
(dotted curve), and 20 mW (dashed curve). 



For example, for 40 mW pump power, the maximum of 
G s and G as are about 1.5 and 0.5, respectively. 



V. SPONTANEOUS GENERATION OF STOKES 
AND ANTI-STOKES PHOTONS: QUANTUM 
CORRELATIONS 

So far we have considered stimulated processes. The 
Stokes and anti-Stokes fields are also generated sponta- 
neously. In this case we have to include input vacuum 
fields. These vacuum fields would be broad band. Thus 
the field at frequency lo s in Fig. [T]is to be replaced by a 
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FIG. 6: (Color online) The normalized output power G s plot- 
ted as a function of the normalized frequency (uj 3 — uj{)/(jO m 
for different pump power, p = 1 mW (solid curve), 6.9 mW 
(dotted curve), and 20 mW (dashed curve). 
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FIG. 7: (Color online) The normalized output power Gas 
plotted as a function of the normalized frequency (u) 3 —uii)/u) m 
for different pump power, p = 1 mW (solid curve), 6.9 mW 
(dotted curve), and 20 mW (dashed curve). 



broad band quantum field Ci n with zero mean value and 
with correlations (6ci n (t)Sc\ n (t')) — 5(t — t'). The cal- 
culations of the output quantum fields are standard [23[ . 
We have used these and introduced the Langevin force 
which stems from the coupling of the movable mirror 
to the thermal environment having zero mean value with 
correlations [2J1 



1 "fm 
7T W m 



2k B T 



j{t ~ t] l + cothf^rr^j dw, 

where ks is the Boltzmann constant and T is the temper- 
ature of the environment. The fluctuations of the output 
fields are obtained as 

6c ovt {u) = V(u)£(u) + E(w)5c m {uj) + F{u)5c\ n {-u), 

(15) 



where 5ci n (u>), and 5c\ n {—uj) are the Fourier trans- 
form of the Langevin force £(t) and the input vacuum 
fields 5cin(t) and 5ct (t), respectively, and where 



F(U>): 

in which 
d(uj) 



2 k 



d{u) 



i[n — i(u> + A)]co, 



{-2^ mX 2 i\cv\ 2 + {u 2 -u 2 m 
1, 



d(u> 

x[k- i(oj + A)]} 

4/ca4x 2 Co ■ 
d(w) ' 

= 4< x 2 A| Co | 2 + (c 2 
x [(k — ioj) 2 + A 2 ]. 



(16) 



(17) 



In Eq. (H3]), the first term containing is the con- 
tribution of the Langevin force acting on the movable 
mirror, while the other two terms come from the input 
vacuum fields. So the fluctuations of the output fields de- 
pend on the Langevin force and the input vacuum fields. 
Further, we define time dependent 5c^ t (t) and 5c^(t), 

where 5c^ t (t) represents the positive- frequency part of 
the fluctuations of the output fields, corresponding to 
Stokes component, and 



6c. 



1 f c 
2t Jo 



Sc out (uj)e lu,t dw, 



(18) 



whereas (t) represents the negative-frequency part 

of the fluctuations of the output fields, corresponding to 
anti- Stokes component, and 



1 

2^ 



Sc out (Lo)e lult cLj. 



(19) 



In the context of Stokes and anti-Stokes radiation gen- 
erated by single atoms, several authors (20l. l25l - |27j found 
important quantum correlations between the Stokes and 
anti-Stokes radiation. Such conclusions were drawn from 
the structure of photon-photon correlations. Motivated 
by these studies and the fact that we are dealing with 
a macroscopic system like a nanomechanical mirror; we 
examine photon-photon correlations in the generated ra- 
diation. 

In the following, like in the work of Kolchin et al. 
[2fj| , we do not differentiate between the Stokes and anti- 
Stokes photons. We calculate the coincidence probability 
defined by 
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9 {2) {r) = 



(0\Scj ut (t)5cj ut (t + r)Sc out (t + r)6c out (t)\0) 
(0\6cl ut (t)6c out (t)\0)(0\5cl ut (t + T)5c out {t + r)|0) ' 



(20) 



in which r is a time delay, and 

" + OO 



1 /■ 

8c out {t) = — / <5c out (cj)e ^du. 



(21) 
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FIG. 8: (Color online) The normalized second-order correla- 
tion function g^'(r) as a function of the time delay t(/is) for 
different pump powers at T = OK. p=l mW (solid curve), 
and 4 mW(dotted curve). 



Now we would evaluate the photon-photon correlations 
of the output fields numerically. We choose the pump 
power p=l and 4 mW and the temperature of the en- 
vironment T = K; the other parameters are the same 
as those mentioned in Sec. III. The correlation func- 
tion g( 2 '(r) between the spontaneously generated pho- 
tons versus the time delay r for different pump pow- 
ers at a temperature of T = OK is displayed in Fig. [5] 
We find that g^ 2 \r) is symmetric. It is also seen that 
g^ 2 \r) > g^(0) as t ^ 0. This demonstrates the pres- 
ence of photon antibunching, which is definitely of quan- 
tum origin. Further, we note the Cauchy-Schwartz in- 
equality g( 2 \r) < g( 2 \0) is violated, and the degree 
of the violation of the Cauchy-Schwartz inequality be- 
comes smaller with increasing pump power. For pump 



power p = 1 mW, the peak value of g^ 2 \r) is about 17, 
and 5 (2) (0) w 3; thus, g (2) {t) / g {2 \0) «5.6. However, 
for p = 4 mW, the peak value of g^ 2 \r) is about 11.5, 
and gW(0) « 3, so g {2) {t) / g {2 \0) «3.8. Therefore, the 
spontaneously generated photons from the optomechan- 
ical system at T — K are correlated nonclassically, and 
the nonclassical correlation becomes weaker with increas- 
ing pump power. This is reminiscent of the parametric 
downconversion process which at low pumping powers 
produces significant quantum correlations. 

VI. CONCLUSIONS 



We have shown that an optomechanical system driven 
by a pump field and a Stokes field can lead to generation 
of a four-wave-mixing signal. The Stokes field is ampli- 
fied. We also find that normal-mode splitting occurs in 
both the generated fields, that is, in both Stokes and anti- 
Stokes fields. We also report lifetime splitting for pump 
power less than a critical power. Further, we have dis- 
cussed the correlations of the photons generated from an 
optomechanical system by spontaneous processes. We 
find the correlations between these photons manifest 
the antibunching effect, and violate Cauchy-Schwartz in- 
equality. Further, the violation of the Cauchy-Schwartz 
inequality becomes weaker with increasing pump power. 
Hence, the optomechanical system can be used to gen- 
erate pairs of photons with quantum correlations. Thus 
the study of both stimulated and spontaneous Stokes and 
anti- Stokes signals provides us with a useful technique for 
studying the strong coupling regime of cavity optomc- 
chanics, as well as quantum fluctuations at macroscopic 
level. 
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